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Abstract 

The spectra which occur in numerical density-matrix renormalization group 
(DMRG) calculations for quantum chains can be obtained analytically for integrable 
models via corner transfer matrices. This is shown in detail for the transverse Ising 
chain and the uniaxial XXZ Heisenberg model and explains in particular their ex- 
ponential character in these cases. 
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1 Introduction 



The density-matrix renormalization group (DMRG) introduced by White in 1992 U is 
a numerical procedure by which one- dimensional quantum systems can be treated with 
spectacular accuracy. It is not difficult to obtain the ground-state energy for spin chains 
with 100 sites up to nine decimal places f2f. Correlation functions can be calculated as 
well and the DMRG has therefore become an important new tool in this area of physics. 

The basis of the method is a proper selection of states in the Hilbert space which are 
important for the target state one wants to study. This is done by dividing the system 
into two parts with corresponding reduced density matrices Q\ and g^. The relevant 
states in each part are then given by those eigenvectors of the g a which have the largest 
eigenvalues. Thus the spectra of the g a enter in an essential way and it is obvious that the 
method will only work well if the eigenvalues drop rapidly enough so that a small number 
of states is sufficient and practically exhausts the sum rules Tr(g a ) = 1. Spectra of such 
a form, where the eigenvalues decrease roughly exponentially, have indeed been observed 
in various calculations [1], [3]. Somewhat surprisingly, however, no detailed investigation 
or discussion of their origin has been given so far. In this article we want to present such 
a study for two particular systems, the Ising chain in a transverse field and the XXZ 
Heisenberg chain. 

Although these systems are integrable and their ground states are known, a direct 
determination of the corresponding g a is difficult. Therefore, we first relate the quantum 
chains to two-dimensional classical systems, namely the Ising model and the six-vertex 
model. As pointed out by Nishino et al. || , the density matrices g a then become partition 
functions of strips with a cut and these can, in turn, be expressed as products of corner 
transfer matrices (CTM's) a la Baxter ||. This feature per se is quite general. For 
integrable models, however, the spectra of such CTM's are known in the thermodynamic 
limit and have, in fact, exponential form, i.e. uj n ~ exp(-an) with integer n. Provided 
the correlation length £ of the system is much smaller than the strip width or chain length 
L, the same should then hold for the g a . This is, indeed, the case: The spectra obtained 
from DMRG calculations are found to agree very well with the CTM results, both in their 
exponential form and in the predicted degeneracies. Deviations occur only at the lower 
end where finite-size and geometry effects are visible. In this way one obtains a simple 
and consistent picture of the density-matrix spectrum and its origin. 

In the following we first describe the relation between density matrix and CTM and re- 
call the results available for the latter. Then we present the numerical calculations for the 
two chains and the comparison with the analytical predictions. The conclusion contains 
some further discussion, also with respect to non-integrable and to critical systems. 
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2 Density matrices and corner transfer matrices 



Consider a spin one-half quantum chain with L sites and Hamiltonian H. The density 
matrix constructed from the ground state |<3>) is 

e =|$)($| (l) 

or in a matrix notation, with cr = {a±, cr 2 , . . . , <Jl] denoting a spin configuration of the 
chain and assuming $ to be real, 

g(cr,cr') = ^(cr)^(cr'). (2) 

By taking partial traces one obtains the reduced density matrices which are of interest 
in the following. Writing cri = {a±, 02, ... , ctm} an d cr 2 = {(Tm+i, &M+2, • • • , one has 

0l(<Tl, O = 51 ®(<7l, 0-2)$(^i, (3) 

and similarly for £> 2 . 

Now imagine that there is a relation of the quantum chain to a two-dimensional lattice 
of classical spins. There are various possibilities for this. For example, H could be a 
derivative of the row-to-row transfer matrix T of the lattice, or T could arise from a 
Trotter decomposition of exp(—/3H). Here we want to assume that the commutation 
relation 

[H,T]=0 (4) 

holds. Then the ground state |$) of H is also an eigenstate of T, and if it also gives 
the maximal eigenvalue (which is the case in our examples), it can be obtained from 
an arbitrary starting vector by applying T a large number of times. Therefore $(<x) 
can be viewed as the partition function (properly normalized) of a half-infinite vertical 
strip with the spin configuration cr at the near end and an arbitrary one at the other, 
far away end. Similarly, g(cr, cr') can be regarded as the partition function for two such 
strips, one extending to —00 and the other one to +00, with end configurations cr and cr', 
respectively. The reduced density matrix Qi, finally, is obtained by identifying <x 2 and cr' 2 
and summing, i.e. by joining the two strips between sites (M + 1) and L, while the rest 
remains unconnected. Thus one arrives at the partition function of an infinite strip with 
a perpendicular cut in it. 

This situation, first discussed in M, is shown in Fig. [I]. The strip in this case is formed 
by a diagonally oriented square lattice. The variables cr\ and cr' x along the cut are shown 
as white circles, the er 2 as black ones. The dashed lines divide the system further and 
define four rectangular corner transfer matrices denoted by A, B, C and D. They are the 
partition functions of the corresponding lattice pieces and of Ramond type, i.e. without 
a common central spin 0. It follows that 
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Figure 1: The strip geometry as discussed in the text with some portions of a square 
lattice. Also indicated are the corner transfer matrices A, B, C, D and the row transfer 
matrix (shaded). 



Ql = ABCD. (5) 

This is the basic relation between the reduced density matrix and the CTM's || Q . 

The CTM's which appear here differ somewhat from the usual ones. In calculations 
for finite-size systems one normally uses the shape shown in Fig. |2] || |TU|, O, |T^, |TB|] . 
The number of spins along both edges is then the same and one is dealing with square 
matrices. This is only a minor point, however, since also the products AB and CD in 
the strip are square matrices for M = L/2. The main point is that, for a large system 
away from criticality, the outer boundary plays a marginal role. In the CTM spectrum it 
affects only the lower end [|K]] and in the thermodynamic limit its effect vanishes so that 



one obtains a well-defined result. 
For this reason the quantity 

q x = ABCD (6) 
where the hat denotes the infinite-lattice limit will have the same spectrum as Qi corre- 
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Figure 2: System composed of four corner transfer matrices with standard shape. The 
arrows indicate the direction of the transfer. 



sponding to the strip, up to some deviations at the lower end and possibly some overall 
shift connected with the different normalizations of the quantities. For the infinite CTM's, 
however, the spectrum is known. Basically, it follows from the Yang-Baxter equations and 
the corner geometry ||. The relevant results will be given in the next section. 



3 Analytical results 

For the treatment of the transverse Ising chain one considers the square lattice of Fig. I] 
with Ising spins at the lattice points and an isotropic coupling K. The transfer matrix 
indicated by shading in the lower part of Fig. [I] commutes with the Hamiltonian []14 , pl| 



# = -E<-<K- A X>X+i (7) 

n=l n=l 

where 5 = cosh 2K and A = sinh 2 2K. The enhancement of the transverse field at the 
right end can be neglected for large L. Then one arrives at the usual homogeneous chain. 

The CTM's for this Ising lattice were studied in and are exponentials of an operator 
which has again transverse Ising form but with site-dependent coefficients which increase 
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linearly with n. Due to the isotropy of the lattice all four matrices are equal and one has 

g 1 = A 4 = e~ Ac ™ (8) 

where 

Hctm = c I £ (2n - 1)< + A £ 2n « +1 1 (9) 

[n>l n>l J 

with a constant c depending on A. The diagonalization of Hctm in terms of fermions 
|16| , [Tl]] then leads to 

oo 

-ffcTM = 51 ( 10 ) 

i=o 



with rij = c]cj = 0, 1 and the single-particle energies 



e 3 



(2j + l)e : A<1 (T>T C ) (11) 
2je : \>l (T<T C ) (12) 



where £ is given by 

e = 'W < 13 > 

Here denotes the complete elliptic integral of the first kind, k' = y/i — k 2 and the 
parameter k with < k < 1 is related to A by 



k 



A : A < 1 (14) 
1/A : A>1. (15) 



An additive constant has been omitted in (|l0|) . From this one finds the following eigen- 
values E for H CTM . 



(a) . A < 1 

The ground state has E = 0. Excited states with one fermion lead to the odd levels 
E = e, 3e, be, . . . With two fermions one obtains in addition all even levels except 2e. 
The level 8e is two-fold degenerate. The three-fermion excitations start at 9e and 
lead to further degeneracies etc. The spacing of the E's, except at 2e, is therefore 
e. 

(b) . A > 1 

The ground state has again E = 0, but it is now twice degenerate due to the long 
range order. This leads to a corresponding two-fold degeneracy of all excited levels. 
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Formally, this follows from the fact that Eq = 0[J The excited states are integer 
multiples of 2e, i.e. the spacing of the _E"s is 2e. 



This is the well-known equidistant spectrum of Hqtm with a spacing which depends 
on A, i.e. on the distance from the critical point. It follows that up to degeneracies, which 
will be discussed later, Qi has a purely exponential spectrum. 

For the Heisenberg model, the associated two-dimensional system is the six-vertex 
model which is specified by three Boltzmann weights a, b, c ||. If the vertex lattice is ori- 
ented parallel to the strip, its row-to-row transfer matrix commutes with the Hamiltonian 



H = Y, R n (16) 

n 

Hn = « + i+*i+A« +1 (17) 

where A is related to the vertex weights via A = (a 2 +b 2 — c 2 ) / (2ab). Strictly speaking, this 
holds only for periodic boundary conditions on T and H. As before, however, the effect 
of the boundaries should not matter for a large system. The non-critical case corresponds 
to A > 1, so that H describes a uniaxial antiferromagnet. 

The infinite CTM's for the six- vertex model are also known ||, [18|, [Hj. They are again 
exponentials of an operator similar to (|i~6"D but with site-dependent coefficients. Thus Q\ 
takes the form (|8|) and 

H CTM = cJ2nH n (18) 

n>l 

where the constant c depends on A. Although the problem does not look like a free- 
fermion one, the eigenvalues of Hctm follow once more from (|T0"D and (12), with s now 
given by 

e = arcosh A. (19) 



Therefore the discussion for case |(b)| above can be taken over and one finds again an 
exponential spectrum for Q\. 



4 Numerical calculations 

In order to check these predictions we have carried out standard DMRG calculations for 
chains with up to 100 sites, keeping between 30 and 64 states in each block. The systems 
were built up via the infinite-size algorithm and Q\ refers to one half of the chain. In 
the following, the eigenvalues uj n of qi, ordered according to magnitude, are shown in 

1 The normalization of this state is discussed in (lql. 
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semilogarithmic plots logo; n vs. n. All cases correspond to correlation lengths which do 
not exceed ten lattice spacings and thus are much smaller than L. 

We begin with results for the transverse Ising chain. Fig. |^ shows an example of 
the spectrum in the disordered region (A < 1). One can see that the first seven levels 
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Figure 3: Density-matrix spectrum of a transverse Ising chain with A = 0.8 and 120 
sites, calculated with 30 states. 



are non-degenerate and equidistant apart from a break between the second and the third 
one. This is exactly the pattern found for the £"s of the infinite lattice in Section [3], case 
[a] . The straight line through six of the levels shows clearly the linear behaviour and 
the level spacing (the slope in a plot liau n vs. n) agrees up to 10~ 3 with the theoretical 
value 2.7565 from fll3|). At the lower end, two-fold degeneracies are visible. Beyond that 
the spectrum becomes flat and no further structure can be seen. There, however, one is 
already at 10~ 14 and near the limits of the accuracy of the calculation. For larger values 
of A, where the slope becomes smaller, the situation is more favorable. 

The situation in the ordered phase is shown in Fig. 0] for A = 1.25. One can see a 
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Figure 4: Density-matrix spectrum of a transverse Ising chain with A = 1.25 and 120 
sites, calculated with 30 states. 



clear difference in the structure of the spectrum: There is no break and the multiplicities 
have changed. Their values (2,2,2,4,4,6, . . . ) are just twice the number Pj of partitions 
of the integers j — 0, 1, 2, . . . and thus correspond exactly to what one expects according 
to Section |(b)| . Only the last multiplet is not correct. We attribute this to numerical 
inaccuracies. 

The two-fold degeneracy needs some comments, though. In some calculations for 
longer chains it was lost and one observed two sets of levels (with degeneracies Pj) shifted 
with respect to each other. This can be traced back to a mixing of the ground state and 
the exponentially close first excited state of H in the numerical procedure, since at the 
same time the expectation value (a^) became non-zero. The mechanism can be discussed 
easily in the limit A > 1, where 



1 
72 



(20) 
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with |$ + ) = | TT ' ' ' T) an d |$-) = | II ■ • • I) being the two ferromagnetic states. This 
leads to 

1 



0-1 — - 
yi 2 



<p+)(<p+\ + \(p-)(ip-\ (21) 



where the states \ip±) now refer to part 1 of the chain. Thus Q\ has two degenerate 
eigenvalues 1/2, the rest being zero. The same holds for the density matrix constructed 
from the first excited state 

|$'}^-L[l<I> + > -!<&->]• (22) 

However, a linear combination a|$) + b\<&') with a 2 + b 2 = 1 leads to eigenvalues (a + b) 2 /2 
and (a — b) 2 /2, i.e. to a splitting. By targeting |$) and |<&') in a DMRG calculation and 
forming such linear combinations one can actually shift the two level sequences deliber- 
ately. 

In Fig. the pure spectra without the degeneracies are shown for several A-values. 
One sees that the slope becomes smaller near the critical point (A = 1), as it should, and 
a comparison of the numbers shows again very good agreement with the value 2e for the 
spacing according to Section || Closer to the critical point, however, a larger number of 
states had to be kept. Also here the spectrum becomes flat around 1CT 14 . 

For the XXZ chain, where only the ordered region and the critical phase exist, the 
situation is even better. Here one can work in the subspace S z = which makes the 
calculation more precise with the same number of kept states. Fig. |^ shows the complete 
spectrum for three values of the anisotropy A. In all cases, the first 10 levels can be 
observed and the degeneracies 2 • (1,1,1,2,2,3,4,5,6,8) can be read off clearly, even 
though some very low levels do not fit perfectly. In comparison to Fig. f| the accuracy 
is much better here so that even at 10~ 18 the structures can be seen. This is a very 
convincing demonstration of the theoretically predicted scheme. The perfect two-fold 
degeneracy indicates that there is no admixture of the excited state here. This is in 
line with the observation that (a%) = in these calculations. Finally, in Fig. [7| pure 
spectra are shown together with linear fits for the first 8 levels. One should note that 
there is an arbitrariness in the definition of degeneracy. Here the criterion was a relative 
difference of less than 10~ 4 . If one relaxes this, the linear region becomes even larger. The 
e-values given in the Figure agree with (|TJ) up to 3-4 decimal places. Also here there is a 
flattening of the spectrum at the lower end. This, incidentally is in contrast to results for 
usual finite-size CTM's, where the spectrum becomes steeper at the lower end |TIJ [TT], [13|1 , 



and can have two origins: The different geometries of the systems in the two calculations 
(Figs. [I] and 0) and limitations due to the truncation in the DMRG procedure. 
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Figure 5: Density- matrix spectra of transverse Ising chains for three values of A without 
the degeneracies, calculated with 30 states for A = 2 and 1.25 and with 60 states for 
A = 1.11. The e- values were obtained from the slopes. 



5 Conclusion 

We have shown for two integrable, non-critical spin chains how the density matrix reflects 
the properties of the associated corner transfer matrices. The exponential character of 
its spectrum can in these cases be understood as a consequence of the star-triangle or 
Yang-Baxter equations.^ The same will hold for other integrable models and one thereby 
obtains a whole class of examples, where the density-matrix spectrum can be obtained 
analytically. In such cases one could estimate in advance, how many states one has to 
keep in a DMRG calculation to get a certain accuracy. Since even for soluble models 
the correlation functions might have to be calculated numerically, this could be useful in 

2 It is amusing to note that such an exponential form arises also for two coupled harmonic oscillators 
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Figure 6: Density-matrix spectra of Heisenberg chains with 98 sites, calculated with 64 
states for three values of A. 



practice. 

One immediate question concerns, of course, non-integrable systems. The connection 
to CTM's still exists, but there are no general results for their spectra. For the three-state 
Potts model, it was shown that the equidistance of the levels of Hqtm is lost for finite 
temperatures |21] . This is consistent with calculations of the £>i-spectrum, which also show 
irregular level spacing, although the overall behaviour is still roughly exponential KI3] . A 



recent investigation of other non-integrable models even indicates a kind of universality 
for the spectra [23]. The fact that the basic features of the CTM's are still similar was 
also used in direct renormalization calculations for these quantities [Q. 

The second question relates to critical systems. Actually a number of the chains 
treated by DMRG belong to this category. Therefore an understanding of the spectra for 
this case would be quite useful. One does not need integrability here, instead conformal 
invariance may be invoked. However, due to the infinite correlation length the shape of 
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Figure 7: Density-matrix spectra of XXZ chains of 98 sites without the degeneracies, 
calculated with 64 states. The e-values were obtained from the slopes. 



the associated two-dimensional system is more important in this case. The CTM spectra 
have been investigated for the usual CTM geometry as in Fig. 



TT|, O, [13] but one needs 



very large systems in order to see the linear spectrum of Hqtm and the logarithmic size 
dependence e ~ 1 / In L predicted by conformal invariance. The same seems to hold for 
the strip geometry. This problem is presently still under investigation. 
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